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Abstract 

Universal quantum information processing requires single-qubit rotations and two-qubit interac- 
tions as minimal resources. A possible step beyond this minimal scheme is the use of three-qubit 
interactions. We consider such three-qubit interactions and show how they can reduce the time 
required for a quantum state transfer in an XY spin chain. For the experimental implementa- 
tion, we use liquid-state nuclear magnètic resonance (NMR), where three-qubit interactions can be 
implemented by sequences of radio-frequency pulses. 
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I. INTRODUCTION 



Quantum computers are capable of solving some còmput ational problems efíiciently for 
which no efficient classical algorithms are known. Examples include the factorization of large 
numbers jï], searching unsorted databases Q], and simulating quantum systems 0,0. While 
this advantage originates from a different scaling behavior compared to classical computers, 
rather than a higher clock speed, the time required for a single gate operation remains a 
critical issue: Reliable quantum computation becomes possible only if a sufficiently large 
number of gate operations can be completed within the decoherence time of the system. 

An important element of many quantum information processing operations is the transfer 
of a quantum state a\0) + (3\1) from one qubit to another We will refer to this process 
as quantum state transfer (QST). We thus discuss a system that is initially in state \^>)i = 
(a|O)+/3|l))^| , 0)i, where the qubit A is in state a|0)+/3|l) and the other qubits in state 
If we denote the QST operation as T, the state transfer from A to B can be represented 
as T : — > \ty)f = \ip)f(a\0) + e^/?|l))B: The final state corresponds to qubit B in 
state a|0) + e^[3\l) and the other qubits in state \i/j)f. A quantum state transfer must thus 
correctly transfer the amplitudes but not necessarily the phases of the state No condition 
is imposed on the state of the other qubits in the system. 

Currently there are three methods that can implement the QST. The first one is quantum 
teleportation proposed by Bennett et al and has been experiment ally realized in optical 
and liquid-state nuclear magnètic resonance (NMR) systems P, 0] . This method is based 
on quantum entanglement and requires quantum measurements. Classical communication 
is also needed if one wants to determine the phase factor in the final state of qubit B. The 
second method is based on swap operations, where T can be represented as T = USji. Sji 
denotes a SWAP gate that exchanges the states of qubits j and l. To realize Sji, one needs 
external operations to control the qubits other than qubits j and l, such as switching on and 
off the couplings between qubits j or l and the other qubits. 

The third method uses a static spin-network Jf],[ïíJ: the qubits are linearly connected by 
Heisenberg interactions. Qubit A is initialized into state ot\0) + /3|1) and the other qubits 
each into state |0). Under the influence of a suitable static coupling network, the system 
evolves such that qubit B ends up in state a\0) + e^/3\l). Unlike the second method, the 
third method does not require spin couplings to be switched on and off, so that it is one 
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kind of quantum computations with the 'always on' interactions that avoids single-qubit 
operations. Hence it is easy to implement in some solid-state systems [12]. In this article 
we concentrate on the third method. For example, the QST can be implemented in a three- 
spin linear chain with the XY- interactions cr^al + Cy&y + (y\o\ + (Ty&y. 

The initial state is chosen as (a|0) + /5| 1))^]00) by setting spin 1 at the location A into 
state a\0) + /3|1) and the other spins into state 1 00) . Waiting for a period of time tç, = 
one obtains the state |00)(a|0) — (3\1))b, which means that the spin 3 at location B is now in 
state a|0) — Both the initial state and the final state are product state. However the 
middle state can be an entangled states. The relation between quantum entanglement and 



13j. In the three spin chain 



the QST in the spin-network has been well discussed in Ref. 
the maximum transfer distance is 2. If one transfers a state over longer distance, one needs 
to design and generate the couplings between spins in a linear chain, or expand the chain 
into a spin network through introducing the additional spins. The details can be found in 
Refs. |q. 

Like other quantum information processing tasks, QST can be effected with a minimum 
set of gates [IJ], typically consisting of single- qubit rotations and CNOT gates that can 
be implemented through two- qubit interactions Q|. An additional possible resource 

]. Effective three-particle interactions exist in some real 



are three- spin interactions 



physical systems, for example in optical lattices constructed of equilateral triangles ^?|. The 
spin Heisenberg chain with three-spin interactions can exhibit interesting phase transition 
phenomena, such as incommensurate phases [!& chiral phase transitions 1201, or a 

nn 

quantum entanglement phase transition[21l. I22j . 

The three- spin interaction that we consider corresponds to a coupling between next- 



nearest- neighbors controlled by the middle spin [17( . It is a rare resource in some quantum 
systems. In this article, we use the three- spin interactions in the spin XY- chain to increase 
the speed of the QST, and quantitatively describe the advantages obtained by using such a 
resource. 

While nature does not provide three-spin interactions between nuclear spins, they can be 
simulated quite readily in liquid-state NMR j^J. For this purpose, one combines the natural 
two-spin interactions of the type J mn a™a™, where ex™ denotes the z- component of the Pauli 
matrix for spin m, and J mn denotes the coupling constant between spins m and n. In this 
work, we use this approach to generate an effective Hamiltonian with variable three-qubit 
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coupling strength to assess the speed-up of the QST due to three-qubit interactions. 



II. XY SPIN CHAIN WITH THREE-SPIN INTERACTIONS 



A. System and hamiltonian 

To test the speed-up of a state-transfer operation by three-spin interactions, we consider 
a three spin XY chain, which is described by the Hamiltonian 

UxY3 = \o x o x + a y a y + a x a x + a y a y ) + -\o x a z a y - o y a z a x ). \\) 

Here, <7 3 x / y / z {j = 1, 2, 3) are the Pauli matrices and we have set H and the coupling constant 
for the two-spin terms to one. To find an analytical expression for the time evolution of this 
system and determine the conditions for state transfer, we write the Hamiltonian (1) as a 
sum of two commuting parts, H X Y3 = C + D, where 

G =(r x a x + a y (T y + -jO x <y z Oy, V = a y a y + cr x a x - -a y cr z a x . (2) 



B. Propagator and transfer speed 



This decomposition shows directly that this Hamiltonian generates a periòdic time- 



evolution: defining k = \J2 + we find C 2 = D 2 = k 2 I and therefore 
U(t) = 



e -itC e -itD 



» .sin(kt) _ r » ,sin(A;í) _ 
[cos(kt)I - i — p^C] [cos(fcí)/ - i — - — -D]. 



k k 
For times t — nr — nn/k (n integer), the propagator returns to unity, U(t) — I 
The matrix representation of the propagator is 
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where c = cos(kt), and s = sin(Atf). 

The propagator generates a state transfer at times íqst = arcsin / ^ : H A > 0, 



8+2A 2 < 

it effects a transfer from qubit 3 to qubit 1, for negative 3-qubit coupling constant in the 
opposite direction. In both cases, the periodicity of the overall evolution implies that the 
reverse transfer occurs at time t = n/k — íqst- The corresponding propagators are 
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(6) 



for the transfer 3 — > 1. 

Figure H shows the dependence of the QST time on the strength À of the three-qubit 
interaction. The overall cycle time r decreases monotonically when a three-body coupling 
is added to the Hamiltonian. However, for A ^ 0, the state transfer is no longer a simple 
SWAP operation, which exchanges the states of qubits 1 and 3, but the transfer beco mes 
asymmetric, requiring different durations for the two directions. While the overall cycle time 
1 — > 3 — > 1 decreases monotonically with increasing |À|, the slower of the two state transfers 
only gets faster than for À = when |A| > 2.71199. 
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C. State transfer 



To demonstrate the state transfer, we set qubit 1 into a superposition state a|0) + 
with the other two qubits in state 1 00) . Applying the forward state transfer Eq. (5) to this 
state gives 

U(t QSTl ^)(a\0) + /3|1»|00> = |00>(a|0> - (7) 

and similar for the reverse transfer. 

If we write this transfer in density operator notation, it reads 
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This result differs when the second and third qubits are initially in different states: 
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The phase of the superposition in the transferred state contains thus information on the 
state of the other qubits. 
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D. Mixed states and parallel implementation 



These different cases can be implemented in parallel by using a mixed initial state 24 1. 
By choosing a — j3 — l/v2 and adding the four initial states in Eqs. (|51 1TTjl . we obtain 



1 1 
1 1 



^ 1 ^ 

10 

10 

1 



-(al + I 1 )®! 2 ®! 3 . 



In the following, we will ignore the unit operator on qubits that are in a superposition state. 
The state transfer acting on this initial state generates then 



U(t Q ST^ 3 y x I 2 I 3 U(t Q ST^ 



(12) 



For related initial conditions, we find 



*s) = z z y 



(13) 
(14) 



U{t QSTl ^V y I 2 I 3 U{t QSTl 

U{t QSTl ^)all 2 I 3 U{t QSTl ^f = I l I 2 a 3 z . 

Obviously the different phases that we found in the state transfer for the pure initial states 
result in the introduction of correlations when a mixed initial state is used. Only if the 
initial state is not a superposition state (Eq. Hlj) . do we find a state transfer that does not 
entangle the transferred state with the other states. 



III. IMPLEMENTATION IN AN NMR QUANTUM COMPUTER 

The nuclear spin system that we use to implement the stepped-up QST has the natural 
Hamiltonian 

3 __ 

H = -7T £ Vitri + ^Ji 2 ay z + £ J 23 cx 2 z a 3 z (15) 
i=i z z 

where Vi denotes the resonance frequency of spin i. Considering this system as a quantum 
simulator of the Heisenberg spin chain described by the Hamiltonian (1), we generate an 
effective evolution (3) by an appropriate sequence of radio frequency pulses. While it is 
relatively straightforward to generate each of the terms of the Hamiltonian (1), they do not 
commute with each other. A sequential generation of the different terms therefore does not 
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produce the correct overall evolution. Two different approaches allow one to generate such 
an evolution: 



Each term is implemented for a very short duration. In this limit, the corresponding 
propagators are close to the unit operator and the noncommuting terms appear only 



in second order 



251. 



• The evolutions U c (t) = e~ itc and U D (t) = e~ itD are written clS cl product such that 
each factor can be implemented directly. 

For the purpose of this paper, we have chosen the second approach. 



A. Decomposing U(t) 

A suitable decomposition of Uc(t) uses the three operators L^. = Ly = <7^cx^/2, 

and = a\a 2 z a'^ j l c 2. These operators can be viewed as the three components of an angular 
momentum vector L c , because they satisfy the cyclic commutation relations [L%, Ly] = iL^ 
and cycl. In terms of these operators, Uc(t) becomes 

Uc (t) = e-^^^ = e- l ^ t)iSj ^^ +Sj ^ L y +cM?) (16) 

where tan 9 C = and the vector n c = (^fs T^f 2 ) cos @c) gives the direction of the rotation 
axis for Uc(t), as shown in Figure|21 Using angular momentum theory, we rewrite this as 

U c (t) = e-^^e^-^e-^^e'^-^ye^. (17) 



In a completely analogous way, we define L% = Ly = a^ay/2, and L\ 



D — 

J z = 

< 7 y< J 'z a x/'2 as the three components of the angular momentum vector L D . In terms of these 
operators, Up becomes 

U D (t) = e- ÜD = e^ïe^'-ïWe^^^eW'-ïWeW, (18) 

where 6 d = vr — 9 C . 

While the two-spin terms L x and L y in Eqs. ()17j) and (|18|) are relatively easy to implement, 
the three-spin terms L z are less straightforward. We re-write them as 



where 77 is an arbitrary real number. Alternatively, we may transform the propagators e tr]L ^ 
and e tr>L ° as 

and 

and use the decomposition of a\a 2 z a\ into one- and two-qubit operators [3| 

The expressions (I20H22|) are identical to the explicit forms of í)19jl. They use only single- 
qubit operations e í<?i<jQ! and precessions under pairwise couplings, e^ * * , which are easy to 
implement experimentally. 

Without loss of generality, we discuss here only the case À > 0. After some simplifications 
flBQ; Eqs. jnj) and (JIHJ) can be represented as 
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(24) 

respectively, and the three-spin terms are implemented according to Eq. (122(1 . 
B. System and pulse sequence 

For the experimental implementation, we used a sample of Carbon-13 labelled 
trichloroethylene (TCE), dissolved in d-chloroform. Data were taken with a Bruker DRX 
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500 MHz spectrometer. We denote the l· l·l nuclear spin as qubit 2 (H2), the 13 C directly 
connected to X H is denoted as qubit 1 (Cl), and the other 13 C as qubit 3 (C3). The param- 
eters of the system and the NMR spectra are shown in Figures El and EJ The difference of 
frequency between Cl and C3 is A^ 13 = 905.3Hz. The coupling constants are J13 = 103.1Hz, 
J12 = 200.9Hz, and J23 = 9.16Hz. Because of the strongly coupled carbons 28] we describe 
the Hamiltonian of the three-qubit system as 
3 

H = É + \ J v.° X A + + | Jx^alal + u\a\ + a z V 3 ). (25) 

Since we use this system to simulate a linear chain with nearest neighbor and three-body 
interactions, we do not use the coupling between qubits 1 and 3, which represent the end of 
the chain. 

Because our quantum register contains only one proton spin, we can implement the 
rotations e 4 x / y by hard 7r/2 proton pulses, which are selective for qubit H2. We denote 
rotations along the ±x or ±y axis as [±| } 2 x / y - The widths of such pulses are so short that 
they can be considered as ideal rotations. Figures El show the actual pulse sequences that 
we used to implement Uc and Ud- 

Implementing spin-selective operations on the carbon spins turned out to be difficult. 
We minimized experimental errors by replacing selective pulses with non-selective pulses 

rï 

and free precession periods 30J, using, e.g., 

e ±if«r? = e^í^'V^e^^' 3 , (26) 

and 

e ±W = e ±i\<ú z é\< e ^\<&\ (27) 

with m = 1 or 3. The n/2 rotations e 4 x i y act on both carbon spins Cl and C3 and were 
realized by hard n/2 pulses. The z-rotations e z ^ a ™ of individual qubits were implemented 
by the " chemical shift rotation" method of Linden et al. íy| . 

The ± signs in Eqs. (J23|) and (J2"2|) refer to two formally different expressions that represent 
the same overall transformation. Implementing both forms and summing over the result 
turned out to be very useful for suppressing experimental artifacts arising from nonideal 
gate operations. When the operations Uc and U d are concatenated, it is possible to combine 
the last operation of Uc with the first of Ud and realize them as a hard pulse [— f ]^' 3 . 
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C. Experimental transfer of (|0) + \l))/y/2 



As discussed in section lTlBI and shown in Figure^ the transfer from qubit 3 to 1 is always 
speeded up by the three-spin interaction for A > 0. We therefore start with this transfer, 
initializing the system to the state ( 1 000) + |001))/v^- To calculate its time evolution, 
we note that, according to Eq. the state 1 000) is an eigenstate of the Hamiltonian, 
C/(í)|000) = |000). Also from Eq. (gj), we find 

(2kc) 2 - (As) 2 . . 2A;cs-As 2 , . 2s 2 + kXcs . . , , 

U(t)\001) = 1 \ k2 [ ' |001) - i p 1010) - 1100). (28) 

We monitor the progress of the state transfer by the amplitudes of the states 1 001) and 1 100) : 
In a superposition with state 1 000) , they correspond to x- magnetization of the qubits C3 
and Cl, respectively. 

As discussed in sections III Cl and III D| we can observe the transfer from the 4 initial 
states |00)(|0) + \l))/y/2, |01)(|0) + \l))/y/2, |10)(|0) + \l))/y/2 and |11>(|0) + \l))/y/2 to 
their respective final states in parallel by preparing their sum as a mixed state I l I 2 a^. using 
pulse sequence 

[^] 2 - [grad} z - [|£ - [grad} z - 



where [grad] z denotes a gradient pulse along z- axis. As usual |32|, we describe these mixed 
states in an operator notation that refers only to the traceless part of the density operator. 
Since [D, a^} = [C, D] = 0, the evolution of this initial condition is determined by C alone, 

Pl (t) = U(t)o*U{t)* = Uo{t)o»U (t)* 

(2kc) 2 - (As) 2 3 2kcs-Xs 2 2 3 2s 2 + kXcs 1 2 3 
= Ik 2 ° x k 2 ay ° z k 2 °* a * a *' ^29) 

The first and last term in Eqn. (j29j) correspond to directly observable magnetization. We 
can therefore monitor the progress of the quantum state transfer by simply recording the 
free induction decay (FID) signal and calculating its Fourier transform. Figures |ü] show the 
corresponding 13 C NMR spectra observed before and after the QST, using TCE. The initial 
condition shows that the signal is concentrated on qubit C3 shown as Figures El (a-c). After 
the transfer C3 — > Cl, the system is in state —ala 2 a z . The main signal is on qubit Cl, 
shown as Figures El (d-f) corresponding to A = 0, 1.5, and 4, respectively. The different 
resonance lines indicate that the magnetization on qubit Cl is aligned along the positive or 
negative x- axis, depending on the state of qubits H2 and C3. This agrees well with the 
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prediction of Eqs. (|SHTT|). After the transfer C3 — > Cl — > C3, the system is in state a^. The 
main signal returns to qubit C3, shown as Figures El(g-i) corresponding to À = 0, 1.5, and 4, 
respectively. The time of QST C3 — > Cl is measured to be t = 1.00, t = 0.62, and t = 0.50, 
and the time of QST C3 — > Cl — ► C3 is measured to be t = 2.00, t = 1.75, and t = 1.13, 
when A = 0, 1.5, and 4, respectively. Here we also use to as the time unit. Compared 
with the case of À = 0, the speed of QST is increased by the three- spin interactions. The 
experimental errors mainly result from the strong coupling between the two carbons and the 
effects of decoherence. Moreover the imperfection of the pulses, especially the 7r pulses for 
refocusing is another error source. 

Similarly the process of transferring (|0) + i\l))/y/2 can be observed by choosing the 
initial state as -f 1 / 2 ^ where a y = (|0) +i|l))((0| — z(l|) — I. Using [C, Oy] = we then have 

P2(t) = U{t)ap{tf = U D {t)ap D {tf 

(2kc) 2 - (As) 2 3 2kcs - As 2 2 3 2s 2 + k\cs 1 2 3 
= 4F °* + P ax ° z P a v a z°z- (30) 

For these initial conditions, it is thus sufficient to consider only part of the evolution 
operator, generating either Uc(t) or Uc{t). 

D. General initial conditions 

For other initial conditions, the full evolution operator U (t) is required. As an example, 
we choose Z 1 / 2 ^ as the initial state, and obtain 

Ps(t) = U{t)alU\t) 

(2kc) 2 - (As) 2 A2kc) 2 - (As) 2 3 2£;cs-As 2 2 3 2s 2 + A;Acs 1 2 31 

= — ^ — ( — — ^ p — w + — p — W^J 

2A;cs - As 2 (2A;c) 2 - (As) 2 2 3 2A;cs - As 2 2 2s 2 + A;Acs l 2 

+ — — i ^ W + — — *z + — p — W 

2s 2 + k\csA2kc) 2 - (As) 2 , „ „ 2A;cs - As 2 , „ 2s 2 + A;Acs 
+ p P \ k 2 vltt p °X + k 2 ( 31 ) 

Noting that J23 is much smaller than J 12 , one finds that Up requires a longer time to 
complete than Uc- For example, when A = 1.5, Uq requires about 340 ms for QST from 
C3 to Cl, while Ud requires about 420 ms. The effective T2 (T 2 ) of the current sample is 
measured to be 0.35s, 0.26s, and 0.23s for Cl, H2, and C3, respectively. When Ud or the 
full U is applied, decoherence results in a significant degradation of the experimental data. 
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We therefore show here only the results of the simulation. For this purpose, we also neglect 
the small strong-coupling effects between qubits Cl and C3. The initial states are chosen as 
J 1 / 2 ^, I l· I 2 ay and l x l 2 a\ as respectively. Because the relevant terms in Eq. (|31|). are not 
directly observable, we apply readout pulses and [|]^ to p%(t) to obtain the observable 
signals of Cl and C3, respectively 

Figures [71121 show the progress of the QST. For each initial state, the results for À = 0, 
1.5, and 4 are given. The data points can be well fitted by the corresponding theoretical 
graphs. Points A, B, and C denote the màxima corresponding to the time of QST C3 — > 
Cl; points D, E, and F denote the màxima corresponding to the time of QST C3 — > Cl — > 
C3. Obviously the time required for the QST decreases with the increase of À. 



IV. DISCUSSION 



s 
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5*1 2 5*23 5*1 2 



(32) 



The QST can also be implemented by a series of SWAP operations. For the three- spin 
chain, the state of spin 1 can be transferred to spin 3 through 

^ 1 0^ 
1 
1 
1 
1 
1 
1 
1 

When À = Eqs.© and (j^J) are equivalent to S13 (up to some phase factors), just as 
discussed in Ref. |26j]. When A ^ 0, however, neither Eq.(jüJ) nor (jSJ) is equivalent to S13. 
The difference between the stepped-up QST and the SWAP operation comes from the three- 
spin interaction, which breaks the symmetry for exchanging spins 1 and 3. One can prové 
that when spin 1 and spin 3 are exchanged, the three- spin terms in Eq. (0) are changed 



V 



from 



A/^1^2,,.3 



nlrr-a:: - a},a 2 al) to -^{cr^a y 



A/^.l_2^-3 



a x a 2 a 3 ) 

y z x) 



Such asymmetry can also explain why 



2\x^z~y "y"z~xl 

ÍQ5T 3 _i differs from t Q sT^ 3 - 

While we have considered here only single-qubit states, it is also possible to transfer 
multi-qubit states through the Heisenberg XY spin chain 0, [jíj, even entangled ones. Such 
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transfers can also be speeded up by three-spin interactions. For example, the four Bell-states 
( 1 00) ± 1 1 1 ) ) / a/2 , ( 1 1 ) ± |10))/\/2 can be transferred from spins 2 and 3 to spins 1 and 2 by 

C/(í Q 5T^ 3 )(|0)i|00) ± |ll)) 23 /v / 2 = (|00) =F |ll)) 12 |0) 3 /v / 2 (33) 

tf(WJ(|0>i|01) ± |10)) 23 /V2 = (|01) ± |10)) 12 |0) 3 (^ + ^^)/v / 2· (34) 

Using the analysis of section |Hj one finds that when À < 0, the speed of transferring the 
entangled states is increased by the three-spin interactions. 

V. CONCLUSION 

We simulated a spin XY chain with three-spin interactions, using a three qubit NMR 
system. Compared to the case where the system contains only two-spin interactions, the 
three-spin interaction increases the speed of the operation. Our results [Eqs. (jüj) and (JSJ)] 
show that when the three- spin interactions exist, the QST is not equivalent to the SWAP 
operation any more. Unlike the SWAP operation, not all rows in the unitary evolution to 
realize the QST have only one nonzero terms. 

The simulation of the XY chain with three- spin interactions offers a possible laboratory 
to study the problems related to three- spin interactions. Our techniques can simulate 
the chain with arbitrary À. In fact À represents the ratio of the three-body and two- body 
coupling constants, because we have set the two- body coupling constants to 1. In a practical 
sense À can be enhanced through increasing the three-body couplings or decreasing the two- 
body couplings to speed up the QST. Although our results are obtained using three spin 
system, they are helpful for the case of more than three spins. 
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FIG. 2: The frame for operations Uc{t) and Uo{t). The vectors n c and n<j denote the directions of 
the rotation axes for the two operations, respectively. They are tilted from the z axis by the angles 
9 C and 6 d = vr — 9 C . The projections of n c and rid into the xy- plane are identical and indicated by 
a black line. The angles between the projection and the x and y axes are tt/4. 
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FIG. 3: The parameters of Carbon-13 labeled trichloroethylene (TCE). The diagonal terms in the 
table are the shifts (in Hz) of the carbons and protons with respect to the reference freqüències 
500.13MHz and 125.76MHz, respectively. The non-diagonal terms are the coupling constants, also 
in Hz. 
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FIG. 4: The carbon spectrum (a) and proton spectrum (b) obtained by applying selective readout 
pulses to the system in its thermal equilibrium state. Each qubit gives rise to four resonance 
lines, which correspond to specific states of the other qubits. The highest frequency lines always 
correspond to the other qubits being in the |00) state, the lowest frequency lines to the |11) state. 



19 



nin 


Ul l U lli 


1 Ulli 


LUi 


Y 1 

n li li i i in 

mi 1 1 Ml 


1' 1 1 Ull 1 I 

1 I I Ull 1 


1 Mil l 

i mi i 


I ï 1 ' ! ! Y 
i mi n 


, 

x 

(d^djd^d,! 1 l 1 1 1 id^djd.l 

; ! ! ! ! kl d =kk! ! ! ! ! 

MIMI 


d, d 4 ldJd 4 Id 4 i i i i 

! kkkk! ! ! ! ! 1 kkk 

II I li 1 lli 1 


i 1 1 1 1 1 1 1 1 
i 1 1 1 1 1 1 1 1 

1 «d Id.idjd,, 1 | 1 

d J 1 1 ! 1 1 kkk 

mill 


! ' 5 ! ' 1 X 

i id^idJd,! 

*1 !! ! ! ! 

III 


Xy y y 

Y 

« stepl 


>X >Y YX 
X XY YX 

-*«step2^« step3^-« step< 


XX Y 
» Y 

t^-« step 


Y X 

5 



b 



C1 


1 1 II II II INI 


II MIM 


\l 


X 
X 

C3 


llll lllíl 


II III 


II II 


im 


llll 1 


í 


i 


1 X 
1 X 


Y 

Y 

H2 


,iiiiiiiiiiii 

, d n l d^d^d,! 1 1 1 1 1 ld 1 ld | |d 1 

! ! ! ! ! W d »J ! ! ! 

MIMÍ 


llll 
llll 

dj 1 1 1 

! kkk 

iiii 

iiii 

Nil 


im i i i ÉT 


iiii 
iiii 
iiii 

d s kkW 

li in 


djdjdjdj i i i 

llll kkk 

III 


d 2 | 

í 


llll 
llll 

llll 
llll 
llll 
llll 

lili 


1 Y 

d i' 

I 


X 


Y Y Y 
Y 

step6 


xx 

x 

- - step 


XY V> 

X Y Y> 

r^»- step8 -- 


5 Y 

■ step9 - - 


Y 

stepl 0- 


Y 1 1 1 


'x 



FIG. 5: Pulse sequences for the implementation of Uc(t) (a) and Unit) (b). Steps 1-10 correspond 
to the ten compound operations separated by 'x' in Uo(t) and Ua(t), respectively. The unfilled 
rectangles denote n/2 pulses, and the filled rectangles denote n pulses. X, X, Y and Y below the 
pulses denote the x, —x, y, and — y directions along which the pulses are applied. Those n pulses for 
which directions are not denoted are refocusing pulses. They are applied in pairs in which the two 
pulses take opposite directions to reduce experimental errors. The durations of the pulses applied 
to H2 and the non-selective pulses applied to Cl and C3 are so short that they can be ignored. 
The selective n pulses for Cl and C3, denoted by the green or blue rectangles, are implemented 
as RE-BURP 2j| and Gauss shaped pulses with 6.2649ms and 2.8252ms durations, respectively. 
The delays are d x = d 2 = T5 ^, d 3 = ^[f - arctan(^)], d 4 = ^{^2 + f + 2vr), 
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FIG. 6: The experimental results demonstrating the QST. The initial state is a^; the corresponding 
spectrum is shown in the left hand column. The results of the QST C3^ Cl are shown as Figures 
(d-f), and the results of the cyclic transfer C3 — > Cl — > C3 are shown as Figures (g-i). The three 
rows correspond to increasing three-qubit coupling strength, A = 0, 1.5, and 4. The time required 
for each transfer is shown in the Figures. At t = íqst 3 _i) the three-spin system is in state — a].a 2 z a% 
and at t = íqst 3 _>i^ 3) the system is in state cnjL 
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FIG. 7: (Color online) Progress of the QST, starting from <r^, for different strengths of the three- 
body coupling. The upper part of the figure shows the overlap of the density operator with the 
target state a\.a 2 z a\ as a function of time. The unit ío of the time axis corresponds to the transfer 
time in the absence of the three-body interaction. The data for À = 0, 1.5, and 4 are marked by 
"*", "+", and "x", respectively. The sòlid lines represent the theoretical results, the individual 
points correspond to the simulated data by setting TCE as the weak coupling system without 
decoherence. Points A, B, and C indicate the màxima corresponding to the transfer times C3^ 
Cl and the points D, E, and F to the transfer times C3 — > Cl — > C3. This clearly demonstrates 
the speedup of the transfer by the three-body interaction. 
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FIG. 8: (Color online) Progress of the QST, starting from (jy, for different strengths of the three- 
body coupling. The data for À = 0, 1.5, and 4 are marked by "*", "+", and "x", respectively. The 
graphs are the theoretical results, used to fit the corresponding data. Points A, B, and C indicate 
the màxima corresponding to the transfer times C3 — ► Cl and points D, E, and F to the transfer 
times C3 -► Cl -> C3. 
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FIG. 9: (Color online) Progress of the QST, starting from a^, for different strengths of the three- 
body coupling. The data for À = 0, 1.5, and 4 are marked by "*", "+", and "x", respectively. The 
graphs are the theoretical results, used to fit the corresponding data. Points A, B, and C indicate 
the màxima corresponding to the transfer times C3 — > Cl and points D, E, and F to the transfer 
times C3 -► Cl -> C3. 
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